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We re-evaluate the decays of the Higgs bosons in the minimal supersymmetric 
• standard model (MSSM) where the tree-level CP invariance of the Higgs po- 

tential is explicitly broken by the loop effects of the third-generation squarks 
Oh. with CP-violating soft -breaking Yukawa interactions. This study is based on 

r~| ! the mass matrix of the neutral Higgs bosons that is valid for arbitrary values 

of all the relevant MSSM parameters. It extends the previous work consider- 
ably by including neutral Higgs-boson decays into virtual gauge bosons and 
those into top-squark pairs, by implementing squark-loop contributions to 
the two-gluon decay channel, and by incorporating the decays of the charged 
Higgs boson. The constraints from the electron electric dipole moment on the 
CP phases are also discussed. We find that the branching fractions of both the 
neutral and charged Higgs-boson decays and their total decay widths depend 
strongly on the CP phases of the top (and bottom) squark sectors through 
the loop-induced neutral Higgs boson mixing as well as the direct couplings 
of the neutral Higgs bosons to top squark pairs. 
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I. INTRODUCTION 



The soft CP violating Yukawa interactions in the minimal supersymmetric standard 
model (MSSM) cause the CP-even and CP-odd neutral Higgs bosons to mix via loop cor- 
rections [IIH^. Although the mixing is a radiative effect, the induced CP violation in the 
MSSM Higgs sector can be large enough to affect the Higgs phenomenology significantly at 
present and future colhders ||l|,|3|,|5HTl|. 



In the light of the possible large CP-violating mixing, we have studied in Ref. all 
the dominant two-body decay branching fractions of the three neutral Higgs bosons based 
on the mass matrix derived by Pilaftsis and Wagner The mass matrix, however, is 
not applicable for large squark mass splitting. In the present work, we re-evaluate all the 
two-body decay modes of the neutral Higgs bosons with the newly-calculated mass matrix 
1^ which is valid for any values of the soft-breaking parameters. Furthermore, we extend 
the work significantly by including the Higgs-boson decay modes containing virtual gauge 
bosons and those into squark pairs and, also by taking into account the squark-loop contri- 
butions to the gluon-gluon decay modes. In addition, we study the dominant decay modes 
of the charged Higgs boson in the presence of the non-trivial CP-violating mixing. 

This paper is organized as follows. In Sec. |I| we give a brief review of the calculation |^ 
of the loop-induced CP-violating mass matrix of the three neutral Higgs bosons. We also 
consider the constraints by the electron electric dipole moment (EDM) on the space of the 
relevant supersymmetric parameters. In Sec. |T| we discuss the effects of the CP phases on 
the neutral Higgs boson decays. The effects of the CP phases on the charged Higgs boson 
decays are discussed in Sec. |V[ Finally, we summarize our findings in Sec. 0. 



II. CP VIOLATION IN THE MSSM HIGGS SECTOR 



The loop-corrected mass matrix of the neutral Higgs bosons in the MSSM can be calcu- 
lated from the effective potential |]T2|,|T3| 
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with V = (pl + al — (pf — a\, g'^ = {g^ + g''^)/4:, and (pi and a. 
the neutral components of the two Higgs doublets: 



1,2) are the real fields of 
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The parameters g and g' are the SU(2)j;^ and U(l)y gauge couplings, respectively, and Q 
denotes the renormalization scale. All the tree-level parameters of the effective potential (|l]) 
such as ml, m\ and = 1^^121 e*^^^, are the running parameters evaluated at the scale Q. 
The potential (|l]) is then almost independent of Q up to two-loop-order corrections. The 
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super-trace is to be taken over all the bosons and fermions that couple to the Higgs fields. 

The matrix Ai in Eq. (|1|) is the field-dependent mass matrix of all modes that couple 
to the Higgs bosons. The dominant contributions in the MSSM come from third generation 
quarks and squarks because of their large Yukawa couplings. The field-dependent masses 
of the third generation quarks are given by 
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where hh and ht are the bottom and top Yukawa couplings, respectively. The corresponding 
squark mass matrices read: 
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where m~, m~ and m~ are the real soft SUSY-breaking squark-mass parameters, A^ and 
At are the complex soft SUSY-breaking trilinear parameters, and is the complex super- 
symmetric Higgsino mass parameter. 



The mass matrix of the Higgs bosons (at vanishing external momenta) is then given by 
the second derivatives of the potential, evaluated at its minimum point 

(01, 02, ai, 02) = ((0i), (02), (ai), (02)) = (t^cos/5, f sin/5, 0, 0) , (5) 

where v = (\/2 Gp)"^/^ ~ 246 GeV. The massless state = ai cos /5 — 02 sin (3 is the would- 
be-Goldstone mode to be absorbed by the Z boson. We are thus left with a mass-squared 
matrix M'\j for three physical states, a (= ai sin + 02 cos 0i and 02- This matrix is real 
and symmetric, i.e. it has 6 independent entries. The diagonal entry for the pseudoscalar 
component a reads: 
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where rriA is the loop-corrected pseudoscalar mass in the CP invariant theories. The CP- 
violating entries of the mass matrix, which mix a with 0i and 02, are given by 
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where g{x,y) = 2 — [(a; + y)/{x — y)]\og{x/y). The size of these CP-violating entries is 
determined by the re-phasing invariant quantities 

mj — mj mi — mi 

which measure the amount of CP violation in the top and bottom squark-mass matrices. 
In the CP-conserving limit, both A^- and vanish, leading to \m\2\ sin(,^ + 612) = 0. The 
definition of the mass-squared m\ and the dimensionless quantities Xt^b, Rt,b and R[ f,, as 
well as the other CP-preserving entries of the mass matrix squared M'h, can be found in 
Ref. [Q. The real and symmetric matrix Ai^fj can now be diagonalized with an orthogonal 
matrix O; 



(10) 

Our convention for the three mass eigenvalues is mni < 1TLH2 ^ ""^//a- 

The loop-corrected neutral-Higgs-boson sector depends on various parameters from the 
other sectors of the MSSM; m^, tan/5, /i. At, A^, the renormalization scale Q, and the real 
soft-breaking masses, mg, m(j, and m^, as well as on the complex gluino-mass parameter 
M~ through one-loop corrections to the top and bottom quark masses |T^. Noting that 
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the size of the radiative Higgs sector CP violation is determined by the rephasing invariant 
combinations Atfie^^ and Abfie^^, see Eq. (||), we take for our numerical analysis the following 
set of parameters: 

\At\ = \Ab\ = 1 TeV , = 2 TeV , 

^Q,u,5 = l^gl = 0-5 TeV , ^ + Arg(/x) = Arg(M~) = , (11) 



under the constraint; 



$ = Arg(A/ie*«) = Arg(A,/ie*«) . (12) 



We vary the common phase $ as well as m^ and tan (3 in the following numerical studies. 
Our choice of relatively large magnitudes of \At^\ = \Ahfi\ enhances CP-violation effects in 
the MSSM Higgs sector [|. 



Our choice of the set of parameters satisfies the necessary condition to avoid a color and electric- 
charge breaking (CCB) minimum in the direction \Q\ = \U\ = 

But, according to the more general study our relatively large values of \At^b\ and \iJ,\ compared 
to those of M~ ~ ~ could give rise to a dangerous CCB minima in the potential which could be 
deeper than the electro-weak vacuum. Therefore, a detailed study of the CCB minima in the 
presence of the non-trivial CP-violating mixing among the neutral Higgs bosons deserves further 
analysis. 
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Before studying the CP-violation effects on tlie neutral and charged Higgs-boson decays, 
it is worthwhile to examine the present experimental constraints on the above parameter 
values (pAl). The charginos, neutralinos, and squarks of the first two generations are suffi- 
ciently heavy for the parameter set ([TI|) . The lightest squark is the lighter top squark ti, the 
mass of which can be as low as 200 GeV for tan (3 = 4 and $ = 0°. The CP-violating phase 
could weaken the LEP lower limit on the lightest Higgs boson mass significantly ||T7| , p!8| . In 
our analysis we show our results when the lightest Higgs-boson mass is above 70 GeV. 

The one-loop effective couplings of the CP-odd components of the Higgs boson to the 
gauge bosons give rise to the electron and neutron EDM's ||T9| at two-loop level. The effects 
could be significant for large tan/? at large \At\, \Ab\ and so that some region of our 
parameter space is already excluded by the present 2a upper bounds on the electron and 
neutron EDM's : \de\ < 0.5 x 10~^^ecm and \dn\ < 1.12 x 10~^^ecm, respectively [^O)- The 
dark-shaded region in Fig. ^ is excluded by the electron EDM constraint in the ($, m^) 
plane for tan/3 = 4 and 10. The unshaded regions give rriHi < 70 GeV. The EDM con- 
straints are avoided only :n the lightly shaded region. The excluded region becomes larger 
for larger tan f3. Even for tan /5 = 4 as shown in Fig. |l|, some parameter space with small 
rriA and large CP-violating phase can be excluded by the EDM constraint. 



It should be noted, however, that the strong two-loop EDM constraints due to the 
one-loop effective couplings of the CP-odd components of the Higgs bosons to the gauge 
bosons may not be valid if there appears cancellation among different EDM contributions. 
Such cancellations may occur between one- and two-loop contributions, or among two-loop 
contributions themselves [|18[^. In fact, the excluded regions of Fig. |l] disappear if such 
cancellation takes place at 50 % (20 %) level for tan/3 = 4 (tan/3 = 10). In this paper, we 
show our results for the whole parameter space (111]). 



III. NEUTRAL HIGGS BOSON DECAYS 

In this section we discuss the phenomenological consequences of the CP-violating Higgs- 
boson mixing on the total decay widths and decay branching fractions of the neutral Higgs 
bosons. We choose tan/3 = 4 throughout this section, which leads larger CP-violating ef- 
fects than the tan (3 = 10 case . 

The most important decay channels of the Higgs bosons are two-body decays into the 
heaviest fermions and bosons because the Higgs couplings are proportional to the particle 
masses. We refer to Ref. for explicit forms of the two-body decay widths of each Higgs 
boson and the relevant interaction Lagrangians. 

The Higgs boson decays into virtual gauge bosons V^*^V^*^ are also important for the 
Higgs bosons of the intermediate mass region, 110 GeV ^ mn ^ 150 GeV |23-p3[. Because 



of the importance of these 3-body and 4-body decay modes, we extend the previous work 
1^ to include the Higgs-boson decays into virtual gauge boson pairs and those into a lighter 
Higgs boson and a virtual Z boson. The partial decay width of Hi \/(*)\/(*) is given by 
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Y(H \/{*)\/W) = G prnj^^Sy {cos P02i + sin/30, 
p: d^-v^)^ elX^/^{uJi,x,y){X{uJi,x,y) + 12xy) 

where = Fy/mv, uji = mj^Jniy, Sw = 2, 6z = 1, and \{x,y,z) = x"^ + y"^ + — 2xy — 
2yz — 2zx. For rriH. > 2my, — 1)^ + e^] can be approximated by 7r/e5(x — 1) and the 
expression of Ref. |^ is recovered for small e. On the other hand, the partial decay width of 
Hi — > HjZ^*^ is given by 

r^H. ^ H,z^'>) ^ Sr<3^ Z'^-^'' ^ . (14) 

^ Sv^vr h a;f7r[(x - 1)2 + e|] ' ^ ' 

where g^ H z — Ouicos [3 O^j — sin /5 02j) — Oij(cos /3 O^i — sin /5 02j). 

For the parameter set ([T^), the decays of two heavy Higgs bosons into sfermion pairs are 
kinematically allowed and their decay widths are given by 

^<"' - f'f"^ = I5^l4/J'^"'(i- '"/,/'"".• • («) 

where Nc = 3(1) for squarks (sleptons). Finally, the Higgs boson decay width into a gluon 
pair including the sfermion contribution to the decay width is given by 

T{H, ^ gg) = {\S!{m%)\' + \Pf{mlj\') , (16) 

where the dimensionless form factors Sf{s) and -Pf (s) are given by 

'''' 

The couplings 5'^//, S'p//; ^ind gjj and the functions Fgj, Fpf, and Fq are given, for example, 
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Figure |^ shows the dominant partial branching fractions of the lightest Higgs boson de- 
cays with respect to the mass for five different values of the CP phase $ : $ = 180°, 
150°, 120°, 90°, and 60°. The lightest Higgs boson mass is found to be less than 120 GeV 
and, as a result, the main decay channels are bb (upper solid line), r+r~ (dashed line), cc 
(dotted line), and gg (lower solid line). We note that the pattern of the branching fractions 
of the main decay channels. Hi bb and r+r", is almost independent of the CP phase. 
However, for larger rriHi, the branching fraction of the decay Hi W^*'^W^*^ becomes com- 
parable to that of the two-gluon channel for $ = 120° and it can be as large as that of the 
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channel for $ = 90° 0. 



On the other hand, the partial branching fractions for the if2,3 decays are very sensitive 
to the CP phase $ and strongly dependent on their masses 171^2 3 as clearly shown in Figs. |^ 
and ^. The upper solid line in each frame is for 7^2,3 ~>- bb and the lower solid line for 
H2,3 ^9 9- The dashed line is for /72,3 t^t^, the dotted-line for 7/2,3 ^ cc, and two 
dash-dotted lines are for the channels H2 W^*^W^*^ (upper line) and Z^*^Z^*^ (lower line), 
respectively. In addition, the thick solid line shows the branching fractions for if2,3 tt, 
the thick dashed line is for 7^2,3 HiZ^*\ and the thick dotted-line is for 7^2,3 Hi Hi. 
Finally, the thick dash-dotted lines are for the titj decay channels. The decay channels for- 
bidden in the CP invariant theories are marked by filled stars in the legend. These channels 
become significant as the phase $ differs from 180°, being the most dominant channels of 
the heavy Higgs bosons for non-trivial values of 

Relegating the detailed descriptions of the 2-body decay channels of the heavy Higgs 
bosons to Ref. it is worthwhile to summarize a few interesting features of the decay 
modes into virtual vector boson pairs and those into top-squark pairs : 

• The decays H2 — > could become the second dominant channel for M^^ ^ 
2M-\v for a large range of nontrivial $. And the branching fraction of the same decay 
channel for the heaviest Higgs boson, H^, — > W^*^W^*\ is also sizable for the similar 
range of 

• The decay channel H2 titi, forbidden in the CP invariant theories, overwhelms 
all the other decay channels for the most range of $, once the mode is kinematically 
allowed. On the contrary, the branching fraction of the decay channel — >■ titi 
decreases as $ differs from 180°, exhibiting a typical two-state mixing between two 
Heavy Higgs bosons for large m^. 

• For the parameter set (|ll]), the non-diagonal decay channels into — > ^1^2 and t2^i 
dominate almost all the other decay channels of H3 once the channel is kinematically 
allowed when $ differs significantly from 180°. The dominance of these channels can 
be clearly seen near the right end of each frame in Fig. ^ with $ = 90° and 60°. 

The dominance of the top-squark channels in the heavy Higgs-boson decays is mainly due 



to the large values of \At\ and in Eq. (|TT]), which enhance the Higgs-boson couplings to 



the top squarks. For the particular parameter set (pAj), two top squarks mix with each other 
(almost) maximally, leading to a significant suppression of the coupling of CP-even Higgs 
bosons to non-diagonal top squarks. 



^It is worthwhile to note that the decay mode Hi 77 with rriHi < 130 GeV is crucial for the 
detection of the lightest Higgs boson at the CERN Large Hadron Collider (LHC) despite its very 
low branching fraction. A detailed analysis of the dependence of this important decay mode on 
the CP phases is to be reported elsewhere [p5|| . 
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IV. CHARGED HIGGS BOSON DECAYS 



Below the th decay threshold, the decay channel W^^*^Hi, involving the 

lightest Higgs boson, constitutes one of the major decay channels of the charged Higgs boson 
together with the decay channel t^u^-. The fermionic decay channels —>■ th and 

T^Vj. are not affected by the CP-violating neutral Higgs mixing. On the contrary, the cou- 
pling of the H^W^Hi vertices depend strongly on the mixing among neutral Higgs bosons 
such that the phase $ can affect the branching fractions of the charged-Higgs-boson decays. 

The interaction of the charged Higgs boson with a neutral Higgs boson and a W boson in 
the presence of the CP-violating neutral-Higgs-boson mixing is described by the Lagrangian 



sin 13 02i - cos /3 + i Ou) H, d'' H+ I + H.c. 



(18) 



in contrast to which the couplings of the neutral Higgs fields to vector bosons are determined 
by the Lagrangian 



jC-h,vv = gmw{cos(302i + sin/^Osi) Hi 



1 



(19) 



Note that if the HiVV couplings are suppressed, the H^W^Hi couplings are enhanced and 
vice versa because of the orthogonality of the mixing matrix O: + 0|j + Og^ = 1. It was 
recently shown [0,0 that the experimental lower limits from the processes e+e^ — > ZHi 
and e+e" — > vvHi on uih^ and tan/? could be significantly relaxed with the suppression 
of the HiVV couplings due to the CP-violating neutral Higgs-boson mixing. In this case, 
the decay channel W^'^*'>Hi plays a crucial role in confirming the existence of the 

CP-violating neutral-Higgs-boson mixing because this mode must be enhanced due to the 
orthogonality relation. 



The explicit forms of the two- and three-body decay widths of the fermionic modes, 
if+ — » t^*^b,T~^h'T-, can be found in Ref. |2^. So, we present in the present work just the 
expression of the decay width of H+ W+^*^Hi: 
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where g^±^^„^ = sin^Osi - cos/^Ogj + zOn, = mr,/mH±, and -f^ = T^/mH±. For the 
charged Higgs-boson mass mH± ^ + the decay width simplifies in the narrow- 
width approximation to 



r(i7+ ^ W+Hi) 



8V2: 



IT 



(21) 



For our numerical analysis of the charged-Higgs-boson decays, we take tan/3 = 4 and 
the parameter set (pi]) as in the analysis for the neutral Higgs-boson decays. Fig. || (a) 
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shows the branching fraction B{H^ W^^*^Hi) on the plane of mH± and The allowed 
region is divided into 4 parts according to the values of the branching fraction; B > 0.4 
(filled squares), 0.1 < B < 0.4 (filled circles), 0.04 < B < 0.1 (filled triangles), and B < 0.04 
(dots). The partial branching fraction B is enhanced with large CP violating phase such 
that this mode could be significant for larger mH± in the CP non-invariant case in contrast 
to the CP invariant case ($ = 180°). The total decay width of the charged Higgs boson is 
shown in Fig. ^ (b) as a function of mj^±. Each line corresponds to $ = 180°(thick line), 
$ = 80°(dash-dotted line), $ = 70°(dotted line), $ = 60°(dashed line), and $ = 40°(solid 
line). In the region with <l> < 80° and mH± ^ 300 GeV, the total decay width is also 
significantly affected by the CP phase $. We show the branching fractions of the charged 
Higgs boson in Fig. ^ (c) and (d) as a function of mH± for two representative values of $: 
$ = 180° and $ = 70°. The thick solid line is for H+ t(*)6, the solid line for H+ t+u^, 
and the dashed line for W~^''*^Hi. We find that the decay channel if+ W~^^*^Hi 

is strongly enhanced as the phase $ differs from 180°, and the channel becomes the most 
dominant decay mode of the charged Higgs boson for ^ 300 GeV when the values of 
$ are less than 70°. 

To summarize, the charged Higgs boson decays into a W boson and a neutral Higgs 
boson are closely related to the neutral Higgs boson decays into gauge boson pairs due to 
the orthogonality of the mixing matrix O in Eq. (|10|). As a result, the neutral and charged 
Higgs-boson decays, significantly affected by the phase are complementary in confirming 
the loop-induced CP violation in the MSSM Higgs sector. 

V. CONCLUSIONS 

Based on the new calculation of the mass matrix of the neutral Higgs bosons which 
is valid for any values of the relevant SUSY parameters, we have re-evaluated the decays 
of the Higgs bosons in the MSSM with radiative Higgs-sector CP violation. The present 
work extends the previous one |^ by including all the neutral Higgs-boson decay modes 
containing virtual gauge bosons and those into top-squark pairs, implementing squark-loop 
contributions to the two-gluon decay channel. We also study the decays of the charged 
Higgs boson in the presence of the non-trivial CP-violating mixing. In addition, we have 
discussed possible constraints from the electron EDM measurements on the CP phase 

We have found that the branching fractions of the main decay channels of the light- 
est Higgs boson, whose mass is less than 120 GeV for tan/3 = 4, are insensitive to 
the CP phase (s) and the Higgs boson mass. However, the sub-dominant decay channel 
Hi v^*^V^*\ which is sensitive to the CP phase $, can become quite significant as the 
CP phase differs from 180°. On the other hand, the decays of the heavy neutral Higgs 
bosons, which are almost degenerate for > ^mz, depend very strongly on the phase. 
Below the thresholds of the decays into top-quark and top-squark pairs, the main decay 
channels of the heavy Higgs bosons consist of bottom-quark and tau-lepton pairs, (virtual) 
gauge bosons, Higgs and (virtual) Z bosons, and Higgs bosons. The relative importance 
of those decay channels depends significantly on the size of the CP-violating mixing. The 
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decay channels Hi ^ tt , titi , t]t-2 overwhelm all the other decay channels once these chan- 
nels are kinematically allowed. The relative importance of these decay modes also depends 
crucially on the the CP phase(s) and the heavy Higgs-boson masses. Finally, the charged- 
Higgs-boson decays into a (virtual) W boson and a lightest neutral Higgs boson become 
dominant when the neutral Higgs-boson decays into gauge boson pairs are suppressed and 
the charged-Higgs-boson mass is less than 300 GeV. 

To conclude, the radiative Higgs sector CP violation induced from the third-generation 
sfermion sectors could alter the patterns of the neutral and charged Higgs-boson decays 
significantly from those in the CP invariant theories. 
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FIG. 1. The regions of ttia and $ allowed by the experimental electron EDM constraint for 
the parameter set ([lT|) are lightly shaded for both tan/3 = 4 (left) and tan/3 = 10 (right). The 
dark-shaded regions are excluded by the electron EDM constraint while the unshaded regions give 
m^i < 70 GeV. 
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FIG. 2. The partial branching fractions of the lightest Higgs-boson decays with respect to the 
mass niHi for tan /? = 4 and five values values of the CP phase $ ; $ = 180° , 150° , 120° , 90° , and 
60°. The upper solid line in each frame is for Hi bb, the dashed line for Hi t'^t~, the dotted 
line for Hi — > cc, and lower solid line for Hi gg. In addition, the upper dash— dotted line is for 
Hi and the lower one for Hi Z(*^Z(*\ 
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FIG. 3. The partial branching fractions for the H2 decay channels with respect to the mass 
for five values of the CP phase $ for tan /? = 4 as in Fig. ^. The decay channels marked by 
filled stars in the legend are forbidden in the CP invariant theories. 
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FIG. 4. The partial branching fractions for the decay channels with respect to the mass 
niHs /3 = 4 and five values of the CP phase ^ as in Figs. ^ and ^. The thick dash-dotted 

line is for the sum of the H3 — > iiii and — > tit2 + ^2^1 modes. The decay channel — > HiZ 
marked by a filled star in the legend is forbidden in the CP invariant theories. 
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FIG. 5. The branching fractions of the charged Higgs boson and its total decay width for 
tan/? = 4 and the parameter set ([Tl]); (a) the branching fraction B{H^ — > W~^Hi) on the plane of 
m}j± and <I>, (b) the total decay width Tjj+ as a function of mjj± for five values of <I> = 180° (thick 
line), $ = 80° (dash-dotted line), $ = 70° (dotted hue), ^> = 60° (dashed line), and $ = 40° (solid 
line), (c) and (d) the branching fractions of the main charged-Higgs-boson decays as a function of 
mu± for <I> = 180° and ^ = 70°, respectively. 
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